Cosmic infrared background from Population III stars and its 
effect on spectra of high-z gamma-ray bursts. 

A. Kashlinsky^ 



ABSTRACT 

We discuss the contribution of Population III stars to the near-IR (NIR) cos- 
mic infrared background (CIB) and its effect on spectra of high-^r high-energy 
gamma-ray bursts (GRBs) and other sources. It is shown that if Population 
III were massive stars, the claimed NIR CIB excess will be reproduced if only 
~ 4 ± 2% of all baryons went through these stars. Regardless of the precise 
amount of the NIR CIB from them, they likely left enough photons to provide a 
large optical depth for high-energy photons from distant GRBs. Observations of 
such GRBs are expected following the planned launch of NASA's GLAST mis- 
sion. Detecting such damping in the spectra of high-z GRBs will then provide 
important information on the emissions from the Population III epoch and lo- 
cation of this cutoff may serve as an indicator of the GRB's redshift. We also 
point out the difficulties of unambiguously detecting the CIB part originating 
from Population III in spectra of low z blazars. 

Subject headings: cosmology: theory - cosmology: observations - diffuse radiation 
- gamma-rays: bursts - gamma rays: theory 



Introduction 



Zero-metallicity Population III stars (hereafter P3) are thought to have preceded the 
normal metal-enriched stellar populations, but because they would be located at high z they 
are inaccessible to direct observations by current telescopes. If massive, they are expected to 
have left a significant level of diffuse radiation shifted today into IR, and it was suggested that 
the cosmic infrared background (CIB) contains a significant contribution from P3 in near-IR, 
both its mean level and anisotropies (see review by Kashlinsky, 2005 and references therein). 
This has recently received strong support from measurements of CIB anisotropies in deep 
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Spitzer/IRAC images (Kashlinsky et al 2005). If P3 are responsible for even a fraction of 
the claimed NIR CIB they would provide a high comoving density of photons all the way to 
the P3 era. In this Letter we analyze effects of such photons on spectra of high-z high-energy 
gamma-ray bursts (GRBs) and blazars that would be observed with the upcoming GLAST 
LAT instrument to 300 GeV. We show that the entire claimed NIR CIB excess (NIRBE) can 
be explained if only ~4% of the baryons have gone through P3 stars. This would result in 
~ 0.1(1 + 2;)^ photons/cm^ whose present day energy is between 1 and A/im. Such photons 
would provide a large optical depth due to photon-photon absorption for GRBs (and other 
sources) at energies that will be probed with GLAST. Detecting this spectral damping in 
forthcoming GRB observations will provide an important test of the P3 era parameters. 



2. CIB from Population III 

Ay wavelengths ^ 10//m the total flux produced by the observed galaxies matches the 
levels of the CIB within its uncertainties, but in the near-IR (NIR) the claimed levels of the 
CIB are substantially higher than the net fluxes produced by galaxies out to flux limits where 
this contribution saturates (see review by Kashlinsky 2005 for details). Fig. 1 shows the CIB 
excess levels (filled circles) over the net flux from galaxies observed in deep surveys (open 
symbols); the caption discusses the details. The excess is significant at 1/xm ^ A < 4/xm, the 
range we term NIR, and its bolometric flux is (Kashlinsky 2005) : 

nW ^ ^ nW 

i^NIRBE = 29 ± 13^ ; FciBexcess(A;iO;um);iO^ (1) 
m^sr m^sr 

At XtlO^m we evaluated the upper limits shown in Fig. 1 described in the caption. The 
wavelengths ^ 10/xm contribute little, so we adopt the value of Fnirbe for what follows. 

It was suggested that the NIRBE is produced by massive P3 stars at high z (> 10) 
(Santos et al 2002; Salvaterra & Ferrara 2003; Magliochetti et al 2003; Cooray et al 2004; 
Kashlinsky ct al 2004). Significant energy release by P3 is suggested from the recent mea- 
surement of CIB anisotropics in deep exposure Spitzer data (Kashlinsky et al 2005). Be- 
cause P3 stars, if massive, would radiate at the Eddington limit, where L oc M, the total 
flux produced by them is largely model-independent (Rccs 1978; Kashlinsky et al 2004). We 
reproduce briefly the argument from Kashlinsky et al (2004): Each star would produce flux 
where is the luminosity distance. Because for massive stars L oc M, the total comov- 
ing luminosity density from P3 is J n{L)LdL oc f^baryon/*!^, where n(L) is their luminosity 
function and /* is the mass fraction of baryons locked in P3 at any given time. In the flat 
Universe, the volume per unit solid angle subtended by cosmic time dt is dV — c{l + z)d\dt. 
Finally, these stars would radiate at efficiency e (~ 0.007 for hydrogen burning). This then 
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leads to the closed expression for the total bolometric flux from these objects: 
3 /G 1 nW 

Here is the mean mass fraction of baryons locked in P3 stars and = is a suitably 

averaged redshift over their era. The total flux is a product of the maximal luminosity pro- 
duced by any gravity-bound object, c^/G, distributed over the surface of the Hubble radius, 
Rh—cHq^, and the fairly understood dimensionless parameters. Prom WMAP observations 
we adopt Qbaryon^^=0.0224 (Bennett et al 2003) and, since the massive stars are fully con- 
vective, their efficiency is close to that of hydrogen burning (Schaerer 2002), e=0.007. 

Requiring that P3 stars are responsible for the flux given by eq. 1 leads to: 

^ , -. 0.0224 0.007 
/3 = (4.2 ± 1.9) X 10-^^37^ 7^ (3) 

^ 'baryon'^ ^ 

Assuming ~10 this is somewhat less than the t5% value suggested by Madau & Silk 
(2005) and considerably less than the >10% value of Dwek et al (2005). Within the NIRBE 
uncertainty, only ^2% of the baryons had to go through P3. This is not unreasonable consid- 
ering that primordial clouds are not subject to many of the effects inhibiting star formation 
at the present epochs, such as magnetic flelds, turbulent heating etc. The only criterion 
for P3 formation seems to be that primordial clouds turning-around out of the primordial 
density fleld have the virial temperature, Tvir, that can enable efficient formation of and cool- 
ing by molecular hydrogen (Abel 2002; Bromm et al 1999). Assuming spherical collapse of 
gaussian fluctuations and the ACDM model from WMAP observations (Bennett et al 2003) 
the fraction of collapsed haloes at z=10 with Tvir >(400,2000)K is (2.6,5) x 10"^ in good 
agreement with eq. 3 as can be derived from Fig. 2 in Kashlinsky et al (2004). 

Eq. 3 was evaluated from 1 to 4//m, but with signiflcant CIB excess flux outside that 
range, /a would increase. However, at wavelengths ^O.l^a^um the high- 2; emissions would be 
below the Lyman break and would be reprocessed to l^m (Santos et al 2002). If < 10, 
which is unlikely in light of WMAP polarization data (Kogut et al 2003), the rest frame 
Lyman break may be redshifted to ^ 1/xm, but the possible extra CIB excess from <l//m will 
be compensated by /a in eq. 3 decreasing with Z3. At longer wavelengths, the CIB excess 
given by eq. 1 can at most increase /a by ^30%. 

The above estimate is subject to two caveats: First, it assumes the NIR CIB at the levels 
given by eq. 1, which were derived assuming a speciflc set of zodiacal light models. The latter 
may carry large systematic uncertainties, which are not included in the formal uncertainties 
in eq. 3. Our estimate of /a is proportional to the CIB excess summarized in Fig. 1 
and eq. 1 and would change accordingly if these are superseded by future measurements. 
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Secondly, this estimate assumes that P3 were massive stars leading to the effective efficiency, 
e, very close to that of the hydrogen burning, e=0.007. If P3 were less massive than ~ 
50 — IOOMq, the effective e of their energy release would decrease by a factor of a few 
(Siess et al 2002) requiring significantly larger values of for a given -Fnirbe- Independently 
of the CIB considerations, if P3 were less massive than ~ 24OM0, their fraction must be 
very small in order not to overproduce the metallicities of the poorest Population II stars 
(Heger et al 2003), leading to CIB fluxes from P3 significantly lower than eq. 1. 

3. Optical depth to photon-photon absorption at high energies at high z: 
appUcation to forthcoming GRB measurements 

If P3 at early epochs produced even a fraction of the claimed NIRBE, they would 
supply abundant photons at high z. The present-day value of 1^=1 MJy/sr corresponds 
to the comoving number density of photons per logarithmic energy interval, dlnE, of 
n-y = —-r-^ — =0.6 cm~^ and if these photons come from high z then ex (1 + zY. 
These photons also had higher energies in the past, ~(0.1-0.3)(1 + z)eV, providing an abun- 
dance of absorbers for sufficiently energetic photons at high redshifts via 77CIB — e'^e" 
(Akhiezer & Berestetskii 1965 ; Nikishov 1962). Stecker & de Jager (1993) have pioneered 
applications of the 7-7 absorption to constraining the present-day CIB from high-energy 
spectra of low- 2; blazars. Madau & Phinney (1996) and Salamon & Stecker (1998) have 
considered effects by evolving normal galaxy populations on potential future intermediate z 
(^0.5) blazars. However, as shown in Sec. 2, P3 stars are likely to have provided a far more 
abundant source of photons at high z to interact with high-energy gamma-ray photons. 

CRBs are the obvious objects whose high-energy emissions would be damped by the 
absorption from the NIRBE photons. This effect was difficult to detect with EGRET because 
of its low sensitivity. Observations by EGRET have detected only 6 GRB's with one of them 
being a record energy 18 GeV photon (Hurley et al 1994). A successor to the Gamma- Ray 
Observatory, NASA's GLAST is to be launched in 2007. Its Large- Area Telescope (LAT) will 
provide significantly improved sensitivity needed for detecting high- 2; GRBs out to ^^=300 
GeV with better than 10% energy resolution above 100 MeV, and its large field-of-view 
should detect ~ 100-300 CRBs/yr. It is expected that ■-^50 of these would have enough 
high-energy photons to measure spectral indices at S >0.1 CeV with uncertainty better 
than 0.1^. We show in this section that the photons produced by the P3 stars should leave 
a detectable signature by damping the high-energy part of the spectra of high- 2; GRBs. 
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We denote with E and S the present day energies of the CIB and GRB photons re- 
spectively and the primes refer to rest-frame energies, e.g. E' = E{1 + z). The photon- 
photon absorption, being electromagnetic in nature, has cross-section ~ that for the Thomp- 
son scattering, ar, it is a{E',S',x) = ^(7^(1 - P^)[2(3{(3^ - 2) + (3 - /?^) ln(i±|)], where 

P — -y/l — and X — cos 9. The cross-section has a sharp cutoff as /3 — > 1, peaks at 

~ |(7t at /9 ~ 0.7, and is a oc /? for [3^0.6. The mean free path of GRB photons in the 
presence of CIB would be (n^c/s^)"^ ~ 0.8((7T/(j)(lMJy/sr//^)(l + z)'^ Mpc. 

The right vertical axis in Fig. 1 shows the product of gtcHq^ and the comoving photon 
number density for given I^,. The figure also marks the regions defined by the photon- photon 
absorption threshold. Discussion in Sec. 2 makes it plausible that, if the NIRBE originates 
from P3 stars, it must have a sharp truncation corresponding to the redshifted Lyman break, 
or ~ O.lz^^m. Thus we assume that, at least at high z^ there are few photons at wavelengths 
that are shifted today to <l/xm. At longer wavelengths there is no observational evidence 
for CIB excess over that from "ordinary" galaxies containing Population I and II, but at the 
same time only spectra of the sources of very high energy 7-rays at very high z would be 
affected by that range of the CIB. Thus it appears that there exists a narrow wavelength 
window of l//m^A^4//m in which the CIB photons can interact with high-£^, high- 2; GRBs 
and probe the emissions from the, so far putative, P3 era. 

We assume a fiat Universe dominated by the cosmological constant and that the NIRBE 
photons originated from P3 at redshifts higher than that of the GRBs, so that oc (1-1-2;)^. 
In this case, the optical depth due to photon-photon absorbtion is: 

dTGRBiS) p (1 + ^) ... . E' dE' 



V^mil + z) + {1 - + z)-^ J-1 Jj^ '"^^ ''^''^'l + z^ E 

where n^dg is the prcscnt-day photon density corresponding to the observed CIB excess. Fig. 
2 shows the resultant optical depth and contributions to it from different z. We adopted 
-E'max=l-24 eV (1 /xm) and the form of n^cm corresponding to the solid line in Fig. 1. 

The onset of r >1 occurs rapidly at S I mlc^ /{leV){l + z)'"^ ~ 261(1 + z)'^ GeV. 
The GRB spectra at these energies should either be strongly damped or there had to be 
only negligible energy releases from the P3 era. Given the high values of r this would still 
hold even if P3 era produced only a small part of Fnirbe, but assumes that the latter has a 
Lyman cutoff redshifted to ~l/i.m today. If the P3 era extended to lower z, the Lyman cutoff 
would occur at the observer wavelength ~ 0.1z3//m and GRB spectra would be damped at 
proportionately lower £. Thus, the location of this cutoff may also serve as an indicator of 
the GRB's redshift. With the advertised LAT energy resolution of < 10% at £^^1 GeV, one 
could determine GRB redshifts from the damping by the P3 photons to better than ~ 5%. 
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How robust is this result? First, photons due to optical counts of galaxies at much 
later times do not affect much the GRB spectra because r ~ zn^aTRn at low z and they 
contribute n^{z = 0)aTRH ~ a few as Fig. 1 shows. Furthermore, unlike the P3 photons for 
which fly oc (l + z)^, those contributed by galaxies would produce a still smaller contribution 
at earlier times compared to P3. Second, even if the NIRBE is not entirely cosmological. 
Fig. 2 shows that P3 emissions should still lead to a very large optical depth, which scales 
as r oc Fnirbe- Thus P3 would likely be the dominant contributors to the optical depth to 
GRBs at high energies. Finally, the magnitude of r at a given bolometric GIB flux should 
not be sensitive to the GIB excess spectrum because the (already narrow) range of 1-4 /im 
available to dampen GRBs at high z is further decreased by (1 + 2;). 



which, when combined with the sudden onset of the optically thick regime, would lead to an 
identifiable damping by the P3 photons. The damping will affect progressively lower energy 
part of the rest-frame GRB spectra as one moves to higher z. At these energies, the GRB 
emission is likely produced by the inverse Gompton and the synchrotron self Gompton com- 
ponents and is expected to be high (or even dominate) for the typical Lorentz factors involved 
(Piran 2004 and references therein); e.g. Dingus et al (1998) construct an average spectrum 
from four EGRET burst and find the differential photon index of driy/dS oc £;^(i-95±o.25) 
out to ~10GcV. It may be sufficient to use GLAST observations of fairly low z GRBs (say, 
z ~2-4 determined spectroscopically in afterglow observations) to establish the existence of 
the GIB from P3 epochs; if positive then the higher z GRBs can be used to further verify 
that it comes from the P3 epochs and/or calibrate the z determinations from the damping. 



Although spectra of more and more distant blazars are now measured with new instru- 
ments, such as HESS ^, these blazars are still too close for unambiguous detection of the 
P3 emissions. The farthest blazars with known spectra are at z ~0. 13-0. 18 and the spectra 
extend to £ <3 TeV (Dwek et al 2005, Aharonian et al 2005). The right panel in Fig. 2 
shows the optical depth of a blazar at 2;=0.18 produced by the 7-7 absorption 1) due to ob- 
served galaxies, 2) due to NIRBE from DIRBE and IRTS, and 3) omitting the IRTS-based 
point at 1.65 //.m and assuming the lower end of the DIRBE-based results. Even if the entire 
NIRBE is correct and originates from P3, the additional damping from it is small and is 
more pronounced only at S >1 TeV, where measurements and interpretation are difficult. 



Fig. 2 shows that the optical depth from P3 is very high, r 




4. Is Population 3 detectable with low z blEiZcirs? 



^http: / /www. mpi-hd.mpg.de/HESS / 
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This is because P3 contribute to the background over a short range of wavelengths longward 
of ~ 1/xm. If P3 were to contribute only a fraction of the NIRBE, their emissions will not 
contribute appreciably to the observed spectra of z ~ 0.1 — 0.2 blazars, but would be seen 
in the high- 2; GRBs. If GLAST collects a large sample of blazars and other AGNs at z^l, 
Fig. 2 shows that they can also be used to probe the emissions from the P3 era. 

This work was supported by the NSF under Grant No. AST-0406587. I thank David 
Band and Demos Kazanas for useful discussions and comments on the manuscript. 
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Fig. 1. — CIB excess (filled circles) and "ordinary" galaxy (OG) contributions (open symbols) 
vs wavelength. The numbers to compute the CIB excess, i.e. observed CIB flux minus the 
OG contribution, are adopted from Fig. 9 of Kashlinsky (2005 and references therein) and are 
discussed at length there (Sec. 5, Table 5 and beyond). Briefly, the net CIB flux is adopted 
from Cambresy et al (2001) at 1.25/im, from Matsumoto et al (2005) at 1.65 fim, from Gorjian 
et al (2000) and Matsumoto et al (2005) at 2.2 fim, and from Dwek & Arendt (1998) and 
Wright & Reese (2000) at 3.5 /xm and from Matsumoto et al (2005) at 4 /xm. The flux from 
OG is taken from HST counts out to 2.2 /xm (open squares from Madau & Pozzetti 2000) and 
from Spitzer/IRAC counts at 3.6 and 4.5 //m (open diamonds, Fazio et al 2004). At lO/zm 
no CIB excess was observed and the levels of CIB are consistent with the net contribution 
from OG. The upper limits on the CIB excess there are shown where net flux from ordinary 
galaxies is known from SCUBA and ISO measurements. The CIB level at 450 and 850 fim was 
taken from Fixsen et al (1998). At 12 and 24 fim we adopted the lowest u]9])er limits on the net 
CIB flux using 7-ray blazar observations (Stanev & Franceschini 1998; Renault et al 2001). 
They are also largely irrelevant for computations of the GRB photons absorption: vertical 
bars with left-pointing arrows show the range where photon-photon absorption is possible 
for the redshifts and energies indicated. The thick hght-shaded sohd and dashed hues show 
the CIB excess spectrum used in computing the optical depth shown in Fig. 2. 
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Fig. 2. — Left: dr/dz vs z for GRB photon energies shown in the panel. Sohd, dotted, 
dashed, dash-dotted and dash-triple-dotted lines correspond to increasing order in S. Middle: 
The net r vs the GRB photon energy for the GRB redshifts shown in the panel. Solid, 
dotted, short-dashed, dash-dotted, dash-triple-dotted and long-dashed lines correspond to 
increasing order in z. The range of redshifts was chosen to avoid overlap between GRBs 
and P3 era; the latter is assumed to have ended by z—10. Only NIRBE from Fig. 1 is 
assumed in the calculations. This assumption is fairly safe at larger z as this component 
gets progressively larger the ordinary galaxies emissions, but at 2; ~ 1 the latter can still 
contribute (Madau & Phinney 1996). Right: Optical depth to photon-photon absorption for 
a source at z=0.18. Line notations correspond to Fig. 1. Dotted line assumes only ordinary 
galaxies measured in deep counts and that their photons originated at z > 0.18. Thick light- 
shaded lines correspond to the NIRBE: sohd line assumes both the DIRBEl- and IRTS-based 
claims at the central points of the measurements and dashed line assumes the "minimal" 
NIRBE with only the DIRBE-based points (i.e. the 1.65 //m point is omitted) and the GIB 
levels corresponding to the lower end of the error bars in Fig. 1. 



